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Examining the evidence for dust destruction in GRB 980703 
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Abstract. The effects that gamma-ray bursts have on their environments is an important and outstanding issue. Dust destruction 
in particular has long been predicted while observational evidence is difficult to obtain. We examine the evidence for dust 
destruction by GRB 980703, in which various inconsistent measurements of the host galaxy extinction have been made using the 
GRB afterglow emission. We construct a spectral energy distribution from nIR to X-ray to measure the extinction at early times 
and compare this with previous findings. We also construct nIR/optical SEDs at intermediate epochs to examine a previously 
reported decrease in extinction. The extinction is very high for a GRB host galaxy. The earliest extinction measurement is likely 
to be lower than previously estimated, and consistent with most later measurements. In a series of SEDs we do not find any 
evidence of variable extinction. We therefore conclude that there is no clear evidence of dust destruction in this case. 
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1. Introduction 

The search for the predicted effects of gamma-ray bursts 
(GRBs) on their environments is key to understanding the 
properties of their host galaxies, found both nearby and in 
the early Universe. Ionisation by the powerful GRB and af- 
terglow emission has now been seen in the form of opti- 
', cal absorption line variability (e.g. GRB 020813, Dessauges- 
Zavadsky et al., 2006; GRB 060418, Vreeswijk et al., 2007) 
and hinted at from several reports of low significance varia- 
tions in X-ray column density for some GRB s, for exa mple 
GRB 01 1121 dPiro et all Eool GRB 050730 dStarling et all 
20051). GRB 050904 ( Boer et all |2003JCampana et al., 2007; 



Gendre et al., 2007) and GRB 060729 (Grupe et al., 2007) and 



fro m comparison of opti cal and X-ray hydrogen column densi- 
ties dWatsonetallr2007l) . 

Optical extinction should also decrease as we ex- 
pect du st to be destroyed by t he GRB jet out t o 10- 
30 pc J5 raine & Salp e"Q Il979t IWaxman & Drainel [2000; 
Fruchter et al., 2001; Perna & Lazzati, 2002). This should oc- 
cur very shortly after the onset of the GRB, but so far only one 
such GRB had potentially shown this effect, GRB 980703, and 
the decrease in extinction would have occur red from approx 
imately 1 day after the GRB onwards (e.g. Vreeswijk et al. 



1999tlHolland et all 1200 lHStarling et aHl2007l) . The evidence 
of this comes from two sources. Firstly the suggestion by 
Castro-Tirado et al.1 d 19991) that the spectral slope estimated 



from their R- and H-band photometric data from around 1 day 
after the burst indicated a redder value than would normally 
be expected for GRBs, translating into Ay ~ 2.2. This is large 



when compared with the average value of A y — 1.07 found 



from afterglow studies at later times. Secondlv. lVreeswiik et al 
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(1999) compiled spectral energy distributions (SEDs) at four 
epochs covering 2.2-5.2 days after the burst and derive a po- 
tential decrease in the intrinsic optical extinction at the 2-3<x 
level. 

There have been optical extinction measurements of 
GRB 980703 by five groups using a number of different meth- 
ods and instrumentation, and measurements of Ay at differ- 
ent times are far from consistent (Figure [1] Table Q}. There 
is, therefore, a need to return to the apparent variability of 
A v in GRB 980703, which if confirmed would be the first 
observational evidence of this kind. We do this via broad- 
band spectral energy distributions, where possible combining 
an X-ray spectrum with early near-infrared (nIR) and optical 
data, and thereby discuss the observational manifestations of 
dust destruction by GRBs. This issue is particularly impor- 
tant to re-exami ne in the current era of rapid follow-up by 
the Swift satellite dGehrels et alll2004l) . robotic telescopes (e.g. 
Liverpool Telescope and Faulkes Telescopes, 2m in diameter) 
and rapid response mode on larger ground-based telescopes 
such as the 4m William Herschel Telescope and the 8m Very 
Large Telescopes. 

2. Re-evaluating the earliest A v estimate 

2.1. Method 

We use a broadband dataset to create a spectral energy dis- 
tribution at early times when a high value of Ay is reported 
from optical data. We to ok the earliest magnitudes reported in 
Castro-Tirado et alj d 19991) (R and H band) which were taken at 
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Time since burst (days) 


Ay 


Method 


Reference 


0.94 


2.2 


nIR/optical colour 


Castro-Tirado et al. 1999 


1.22 


1.5±0.1 


nIR/optical/X-ray SED continuum 


Vreeswijk et al. 1999 


1.33 


0.9±0.1 


nIR/optical/X-ray SED continuum 


Starling et al. 2007 


3.22 


1.38+0.35 


nIR/optical SED continuum 


Vreeswijk et al. 1999 


4.4 


0.3±0.3 


host spectrum B aimer decrement^ 


Djorgovski et al. 1998 


4.42 


0.84±0.29 


nIR/optical SED continuum 


Vreeswijk et al. 1999 


5.22 


0.90±0.25 


nIR/optical SED continuum 


Vreeswijk et al. 1999 


5.3 


0.9±0.2 


nIR/optical SED continuum 


Bloom et al. 1998 



Table 1. Measurements of optical extinction gathered from the literature. Errors are either given at the \cr level or assumed to be 
so where no confidence level is stated. 
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Fig. 1. Optical extinction (Ay) values reported in the literature 
(squares). All errors are plotted at the \<j confidence level. 



0.94 days since burst (and used for their estimates of A v ). We 
converted the observed magnitudes into fluxes assuming the 
Johnson filters, after applying a correction for the Galactic ex- 
tincti on of E(B - V) = 0.057 dSchleeel. Finkbeiner & Davis . 



due to the redshift of z - 0.9661 (Djor govski et al 



19981) . The Lyman-a Forest correcti on (e.g. Mad au. 995) 



1998; 



Vreeswijk et all 1 19991) . is negligible for these red bands. 

We then calcula ted the required end-time for the BeppoSAX 
(Boe lla et al.L ll997) X-ray observations in order to obtain a log 
observation mid-time of 0.94 days, given that the observations 
began on 0.827 days (Low Energy Concentrator Spectrometer, 
LECS, covering 0.1-4 keV) and 0.844 days (Medium Energy 
Concentrator Spectrometer, MECS, covering 1.4-10keV). We 
extracted a MECS spectrum with these time intervals (up to 
0.96 days) and found that only of order 80 counts were present. 
This did not allow us to create a spectrum of sufficient quality. 



The observed 2- 10 keV X-ray flux a t 1 day is derived from the 
light curve in lGendre & Boerl d2005h to be (4.8+0.07) xl0~ 13 
erg ctrT 2 s -1 with a decay index of ax - 0.9+0.2, (where 



Fy OC f 



1 V -P 



oc t 



-(T-l) 



). We note that the early optical decay 



index of a = 0.85 + 0.84 given in Zeh, Klose & Kannl (120061) 
is consistent with this (though note the large associated error). 
We used this information to transform the full Bepp oSAXLECS 
and MECS X-ray spectra presented in I Starling et alj ( 12007b 
(log observation mid-time 1.33 days) to the time of the early 
optical photometry. 

We fitted the optical and X-ray data together in a coun t 
space fit following the method outlined in Starling et alJd2007l) . 
The underlying continuum is fit with an absorbed power law. 
Previously, wh en taking into account consistency with the 
fireball model (M eszaros & Reesl Il992l I1997J: I Wiiersetal , 



19981: IWiiers & Galamal Il999t iMeszaros & Reesl 119991) . sev- 
eral authors could estimate the position of the cooling break, 
v c . T his break appears to move to lower frequencies with 



time dVreeswijk et all [19991: iBloom et allll998l:IStarling et al . 



2008). It follows from these results that at 0.9 days since 
burst we expect a single power law shaped spectrum spanning 
the nIR/optical to X-ray regime. Galactic abso rption is taken 
into a ccount using Nn = 4.98 x 10 20 ctrT 2 (Kalberla et al 



2005 ). Galactic extinction is fixed at E(B - V) -■ 0.057 
dSchlegel. Finkbeiner & Davisl Il998l) . We modelled the ex- 
tinction in the nIR/optical bands using the three most well- 
known extinction curves of the Milky Way (MW), Large 
Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC), 
as parametrised bv lPeil dl992l) . 



1 We note that this method probes a different volume than the other, 
line-of-sight, methods listed and should therefore not be directly com- 
pared. For a d i scussion on met hods of est imating extinction se e e.g. 
iMathij Jl983l) . ISavaglio&FaHl J2004l) and lStarling etall J2007I) . 
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Model 


r, 


Ay 


xl0 22 cnT 2 


^/dof 


PL+MW 


2.05±0.02 


1 74+0.10 


0.8*}* 


32/24 


PL+LMC 


2.05±0.02 


1.14+0.10 


0-8t1 


32/24 


PL+SMC 


2.04±0.02 


1.08+°; 09 


0.7^ 


33/24 



Table 2. Absorbed power law fits to the SED at 0.9 days since 
burst. All errors are quoted at the lcr confidence level. 



2.2. Results 

The results of SED fits are given in Table [2] A single ab- 
sorbed power law is a statistically acceptable fit to these 
data with power law photon index F ~ 2.0, X-ray absorp- 
tion of order 7xl0 21 cirr 2 and optical extinction of E(B - 
V) ~ 0.4 (approximately equivalent to Ay ~ 1). Comparing 
Table [2] with Table [TJ shows that the extinction values are 
consistent at the lcr level with tho se derived from fits a t 
3.22, 4.42 and 5.22 days since burst d Vreeswiik et all 1 1999ft . 
1.33 days s ince burst ( Starling et al., 2007) and at 5.3 days 
since burst ([Bloom et all 1 1998b . At this confidence level the 
result is inconsistent with previous estimates at 0.94 days 



dCastro-Tiradoetallll999l). 1.22 days (tVreeswijk et all [1 999) 



and 4.4 days (IDjorgovski et all 119981 host galaxy measure- 
ment). MW, LMC and SMC extinction curves are indistin- 
guishable for these data; this is partially because we do not have 
data in the bluer bands where the 2 175 A bump characteristic 
of only the MW curve would occur at this redshift. The SMC 
extinction curve is generally the best approximation to GRB 
host galaxy extinction laws of the three aforementioned curves, 
as shown in many studies and most likely following the low 
metallicity pattern of GRB hosts (e.g . Galama & Wiiers, 2001 
Stratta et all 120041 [Kann. Klose & Zehl 120061: ISchady et al. 



20071 1 Starling etal ll2007h 



3. Re-evaluating variability claims 

3.1. Method 

At four later epochs spanning 1.2 - 5.2 days since burst 
Vreeswiik et al] (1 19991) created nIR/optical SEDs and find a de- 
crease in the intrinsic optical extinction of AAy = 0.24 - 0.96. 
In order to test the possible decreasing extinction proposed by 
these authors we recreated the nIR/optical SEDs at the epochs 

3.2, 4.4 and 5.2 days since trigger. We do not reproduce the 
SED at 1.2 days because this epoch (optical data from 1.2 days 
extrapolated to 1.3 da ys) has been covered using the same SED 
method and models in lStarling et al.l (12007b and we will use the 
results from that fit . Likewise, the previously reported epochs 
of 5.2 and 5.3 days d Vreeswijk et all 1 1 999l Isioom et all 1 1 998 . 
respectively) are close in time so are tr eated here with a single 
SEP. We gathered all ava ilable data (IVreeswiik et all Il999t 
Castro-Tirado et all Il999l and references therein) to create 
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Fig. 2. Unfolded SED (instrumental response at X-ray energies 
removed) at 0.9 days since burst (upper panel) and residuals 
to the fit (lower panel). The black crosses (+) are the optical 
photometry and MECS spectrum and the black stars show the 
LECS spectrum. The solid (red) lines are the model of a power 
law absorbed in X-rays and in optical by an SMC-like extinc- 
tion. 



light curves in seven photometric b ands, BVRIJHK (the ff 
band was not included in the SEDs of iVreeswiik et al.l dl999ft ). 
interpolating from the six nearest points in time to find the mag- 
nitude per epoch. We corrected for the host gala xy contribu 



tion to each band using the host magnitudes from lSokolov et al 



(2001) (B band) and lVreeswiik et al.l (1999) (VRIJHK) and er- 
rors on the host magnitude determinations were combined with 
the datapoint and fit errors. 

3.2. Results 

The SEDs are shown in Fig. [3] and fit results are presented in 
Table [3] Fitting a model consisting of a power law plus fixed 
Galactic extinction and variable (SMC- or MW-like) host ex- 
tinction, we find we cannot constrain both the power law slope 
and the host extinction simultaneously. We derive upper lim- 
its on Ay which are consistent with each other for these three 
epochs, and with most previous measurements. 

We may hope to reproduce the de crease in Ay if we keep 
the power law slope fixed, as done by I Vreeswiik et alj (119991) . 
They fixed B opt = 1.013 (or F = 2.013) which they measured 
from a single power law fit to their optical-X-r ay SED at 1.2 



days since trigger. At a similar time, 1.3 days. IStarling et al. 



(l2007h measured B opt = 0.87 (SMC) and 0.88 (MW) (or 
T = 1.87, 1.88) from a broken power law fit to an optical- 
X-ray SED, so we tested both fixed values for the power law 
slope. Making this assumption, that the power law slope does 
not change in the optical regime over this time period, we con- 
sistently find values of Ay ~ 0.8 - 1.2, in agreement with the 
values reported in the initial study for the final two epochs 
but lo wer than those of the first two epochs dVreeswiik et al. , 
1999). Of those tested, no one extinction model is preferred. 
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Energy (keV) 

Fig. 3. nIR/optical SEDs at 3.2 (upper), 4.2 (middle) and 5.2 
(lower panel) days since burst with best-fitting power law 
model plus SMC extinction overlaid (solid lines). The SEDs 
are plotted in the same units as the optical-X-ray SED of Fig. 
|2]for comparison. 



4. Discussion: is there evidence for dust 
destruction in GRB 980703? 

Construction of a broadband SED at early times (0.9 days) re- 
sults in a lower value for the optical extinction than previously 
found using the optical colour, bringing this first extinction es- 
timate into line with all later afterglow measurements at the 2<x 
level or better (Fig.Q]). We then addressed the possible decrease 
in extinction claimed for several later epochs by creating new 
SEDs at these times with both increased wavelength coverage 
and fitting with a wider range of extinction models. It appears 
that given the substantial errors on most afterglow-derived A y 
measurements they are all consistent when the SEDs are recre- 
ated in a consistent manner and fit with the same models. We 
have tested the need for decreasing extinction across all our 
newly derived Ay measurements spanning 0.9 - 5.2 days since 
trigger by performing a linear fit to the five values and their 
lcr errors. We adopted the SMC extinction curve in all cases, 
with the power law photon index in the optical regime as mea- 
sured at 0.9 and 1.2 days and fixed at 2.013 at 3.2, 4.4 and 5.2 
days (all consistent with F = 2.0). Throughout this time period 
the optical spectral slope should be well described by a single 
power law. The cooling break is proposed to lie in the medium 
energy X-ray band e arly on (e.g. at ~8xl0 17 Hz at 1.3 days 
(Star ling et all 120081) ) to account for the X-ray temporal slope 
within the fireball model (hence a broken power law model was 



Epoch 


Model 


r a,. 


,T/dof 


3.22 


PL+SMC 
PL+MW 
PL+SMC 
PL+MW 
PL+SMC 
PL+MW 


2.6^! <L1 

2.013 0.82!^ 
2.013 0.87^ 

1.87 0.94^ 

1.88 0.99!°;| 


9.4/4 
9.4/4 
9.7/5 
10.2/5 
9.9/5 
10.4/5 


4.42 


PL+SMC 
PL+MW 
PL+SMC 
PL+MW 
PL+SMC 
PL+MW 


2 l*° s <2 2 
3.0+0.3 <2.3 
2.013 0.97^| 
2.013 1.05^ 

1.87 1.08^ 

1.88 1.18+0.28 


8.1/4 
8.7/4 
8.1/5 
9.2/5 
8.1/5 
9.3/5 


5.22 


PL+SMC 
PL+MW 
PL+SMC 
PL+MW 
PL+SMC 
PL+MW 


3.0+0.4 <1.8 

3-0-oi <L9 
2.013 0.88!°;^ 

2.013 0.96!°;*° 

1.87 FOO^I* 

1.88 1-09!°;^ 


1.5/4 
1.55/4 
1.9/5 
2.7/5 
1.99/5 
2.8/5 



Table 3. Absorbed power law fits to the SEDs at epochs 3.2, 4.4 
and 5.2 days since burst. All errors and upper limits are quoted 
at the lcr confidence level. 



adopted for the broadband SED at this time) and in the sof t X- 
ray band or UV five days since trigger (IBloom et all Il998h . In 
a homogeneous circumburst medium the cooling break would 



move to lower frequencies with time approximately as t up 
until the onset of any jet break, and would remain within the 
X-ray band throughout the epochs covered by this study. 

We find that a single value of Ay - 0.96+0.07 is an ac- 
ceptable description of the data, with a fit statistic of ^ 2 /dof 
= 2.13/4. This is shown in Figure |4] A fit to the data using a 
polynomial of degree 1 (allowing for an increase or decrease in 
the extinction) returned a fit statistic of ^ 2 /dof = 1.62/3, which 
does not significantly improve the fit; the F-Test probability for 
inclusion of the additional free parameter is 0.4, i.e. there is a 
40% chance probability of obtaining this result. While the lack 
of broadband (X-ray and UV) data at these epochs prevents 
us from definitively ruling out variable Ay, the results strongly 



R.L.C. Starling: Examining the evidence for dust destruction in GRB 980703 



suggest that a decreasing extinction is not required. This may 
be due to a combination of the additional data used in our study, 
application of the SMC extinction curve and a different fitting 
procedure. We cannot identify the preferred extinction model in 
any of our fits, again because more optical bands, particularly 
towards the UV where the curves deviate most significantly, 
are needed. This also means that where the extinction curves 
applied to estimate the previous Ay values differ, this is un- 
likely to affect the overall results. iKann. Klose & Zehl (|2006) 
have shown that it is difficult to distinguish between optical ex- 
tinction models forz < 1.5 of which GRB 980703 (z = 0.9661) 
is an example. 

Most models predict some form of ionisation and dust de- 
struction by the GRB on its immediate surroundings. To ob- 
serve this we need to probe distances of perhaps u p to 10 - 30 
parsec from the GRB (Wax man & Drainei 120001) . Decreasing 
X-ray column densities, while not yet measured at a highly sig- 
nificant level, suggest ionisation by the GRB is taking place. 
Comparison of H I column densities from optical spectra with 
equivalent hydrogen columns measured from X-ray spectra 
concluded that the optical and X-ray gas are likely not co- 
located, and the wide s pread in measured values implies that 



ionisation plays a role (IWatson et all 120071) . These authors es 



timate a maximum distance out to which most of the hydro- 
gen has been ionized by the GRB of 3 parsec. Time resolved 
optical spectroscopy has now revealed compelling evidence of 
UV pumping, seen as variations in opti cal absorption lin e prop- 
erties, most evident in GRB 060418 dVreeswijk et all I2007J) . 
These results suggest that the optical regime is probing gas that 
lies 1 .7 kiloparsec from the GRB site, further than the distance 
out to which dust is thought to be destroyed. 

Dust destruction will also only take place as the GRB ion- 
ising flux reaches the cloud, perhaps only minutes or seconds 
(in the observer frame) after the GRB has occurred depend- 
ing on the distance to the cloud and the number and energy of 
ionising photons, and should be temporally coincident with the 
ionisation. To return to the example of GRB 060418, we note 
that observations began 1 1 minutes after the burst. Changes in 
the X-ray column density thought to be due to ionisation by the 
GRB are on similar timescales, for example in GRB 050730 
this occurred over the first 8 minutes since trigger. In the case 
of GRB 980703 accurate localisation occurred only >20 hours 
after the burst when an X-ray aftergl ow was found using the 



BeppoSAX narrow field instruments (IGalama et all Il998l) . It 



seems unlikely, therefore, that we would be able to detect any 
observational manifestation of either ionisation or dust destruc- 
tion by the GRB jet in this source. Swift observations of GRBs 
have revealed later time X-ray flaring (obse rved up to 73 hours 
after the onset of the GRB in rare cases ICurran et all 120081) 
which may also have an effect on the GRB surroundings and 
at later times than the GRB jet itself. The X-ray light curve for 
GRB 980703 is, however, not well sampled enough to allow 
flares to be clearly recognised dGendre & Bo er. 2005). 

The lower average optical extinction value we have de- 
rived for GRB 980703 remains the highest of all the extinc- 
tion values found in the BeppoSAX GRB afterglow sample 
dStarhng et all 120071) . In the Swift era even higher extinction 
values have been proposed for so-called 'dark' bursts (e.g. 



I.b 

1.4 
1.2 


[] 






- 


1.0 


- 








- 














0.8 

0.6 

4 


[] 


- 



J 4 b 6 

days since trigger 

Fig. 4. Optical extinction (Ay) values derived in this study from 
afterglow SED fits (squares). The dashed line shows the best fit- 
ting constant value to the distribution and the dotted lines show 
its associated error. All errors are plotted at the lcr confidence 
level. 



Jakobsson et al., 2004; Rol et al., 2005) at similar redshifts to 



GRB 980703, such as GRB 051022 located at z = 0.8 in which 
our line-of-sight to the GRB may pass through a du sty region 



of the host galaxy unrelated to the GRB location dRol et al 
20071). 



5. Conclusions 

The optical extinction measurements in the literature for 
GRB 980703 are inconsistent, and suggested a decrease in ex- 
tinction with time indicative of dust destruction which is pre- 
dicted to occur in GRBs. We have re-examined the greatest out- 
lier in these measurements, using broadband SED fits rather 
than optical colour, to estimate the extinction. We find a lower 
value than previously estimated, in line with most of the later 
measurements. We have also investigated the apparent decrease 
in extinction by recreating the nIR/optical SEDs at these epochs 
with extended wavelength coverage and a wider range of mod- 
els for the extinction law. Both from our new SED fit results and 
given that the extinction measurements begin as late as 0.9 days 
after the burst, we conclude that there is no evidence of dust de- 
struction by this GRB. In the current era of rapid accurate GRB 
localisations and follow-up we now have an opportunity to ob- 
serve some of the effects of GRBs on their environments, and 
studies such as this using multiple early-time SEDs will be of 
great value. 
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